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ABSTRACT 

X-ray data of the Centaurus cluster, obtained with XMM-Newton for 45 ksec, 
were analyzed. Deprojected EPIC spectra from concentric thin shell regions 
were reproduced equally well by a single-phase plasma emission model, or by 
a two-phase model developed by ASCA, both incorporating cool (1.7-2.0 keV) 
and hot (~ 4 keV) plasma temperatures. However, EPIC spectra with higher 
statistics, accumulated over 3-dimentional thick shell regions, were reproduced 
better by the two-phase model than by the singe-phase one. Therefore, hot and 
cool plasma phases are inferred to co-exist in the cluster core region within ~ 70 
kpc. The iron and silicon abundances of the plasma were reconfirmed to increase 
significantly towards the center, while that of oxygen was consistent with being 
radially constant. The implied non-solar abundance ratios explains away the 
previously reported excess X-ray absorption from the central region. Although 
an additional cool (~ 0.7 keV) emission was detected within ~ 20 kpc of the 
center, the RGS data gave tight upper limits on any emission with a tempeartures 
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below ~ 0.5 keV. These results are compiled into a magnetosphere model, which 
interprets the cool phase as confined within closed magnetic loops anchored to 
the cD galaxy. When combined with so-called Rosner- Tucker- Vaiana mechanism 
which applies to solar coronae, this model can potentially explain basic properties 
of the cool phase, including its temperature and thermal stability. 

Subject headings: conduction — magnetic fields — plasmas — galaxies: clusters: 
individual (The Centaurus Cluster) — X-rays:galaxies:clusters 



INTRODUCTION 



Intra-cluster medium (ICM), i.e., the hot plasma confined within the gravitational po- 
tential of clusters of galaxies, constitutes the most dominant known form of baryons. The 
ICM was thought to cool over the Hubble time by emitting optically-thin thermal X-rays, 
because its radiative cooling time is estimated to be ~ 10 8 yr at the center of many "cD 
clusters" (those hosting cD ga laxies at the ir centers). This inspired so-called cooling flow 
(hereafter CF) hypothesis (e.g. lFabianlll994j ). The idea was apparently supported by several 
soft X-ray features observed around cD galaxies, including general decreases of the ICM 
temperature, strong excess X-ray surface brightness, and excess photoelectric absorption up 
to a few times 10 21 cm -2 . 

With the first imaging capability over a broad X-ray energy ba nd up to 10 keV and a 



much improved energy resolution than was available before, ASCA ( ITanaka. Inoue fc Holt 



19941 ) has provided a number of new results that altogether cast serious doubt upon the 
reali ty of CFs. Thes e include; a shortage of the predicted cool gas in the Hydra-A clus- 
ter (jlkebe et al.lll9970; the presence of uncooled ICM in 3 -dimensional core regions of the 
Centaurus cluster (IFukazawa et al.lll994l ; llkebe et al.l 119991 1: a central excess brightness of 

hierarchical 



Abell 1795 that does not depend strongly on the X-ray energy 



potential structures in several galaxy groups (jlkebe et al .111996 



Xu et al. 



1998 



Matsushita et al.lll998l ); and 



systematic differenc es in the ICM chemic al composition between the cent ral and outer re 



gions of cD clusters (IFukazawa et al.ll2000l ). We compiled all these results in lMakishima et al. 



( 120011 . hereafter Paper I), and argued that the classical CF hypothesis needs a significant 
revision. The AS CA suggestions have been reconfirmed and reinforced by Cha ndra and 
XMM-Newton (e.g. iTamura et aDboOlal : iKaastra et aDboOll : IPeterson et aDboOlh . 



Although the CF scenario in its original form is no longer considered valid, the ICM 
temperature of a cD cluster does decrease toward the center. Furthermore, the deprojected 
radial ICM temperature profiles of a fair number of clusters all reduce to a single "univer- 
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sal" profile (jAUen. Schmidt fc Fabian!l200ll ; IKaastra et al.ll2004l ). which reach, at the cluster 
center, a minimum value of 

T c = (1/2 ~ l/3)T h (1) 

jAUen. Schmidt fc FabianlhoOll : IKaastra et all hooj ). Tere, Th is the ICM temperature ob- 
served outside ~ 100 kpc. It is yet to be explained how the intriguing universal temperature 
profile and the scaling of eq.([T]) are produced, and how they are related to mechanisms which 
suppress CFs. 

The concept of the universal temperature profile is based on a "single-phase" (hereafter 
IP) picture, that the ICM at a given three-dimensional radius is represented by a single 
dominant temperature which monotonically decreases toward the center. In contrast, the 
original ASCA results on cD clusters (Paper I) are based primarily on a "two-phase" (here- 
after 2P) view; the ICM consists of a "hot phase" and a "cool phase" , intermixed together, 
with the volume filling factor of the latter increasing toward the center. This view was at 
the beginning motivated by the complex angular response of the ASCA telescopes, but was 
reinforced a posteriori by the fact that 2P fits to the ASCA spectra of a cluster generally give 
a pair of well-defined temperatures regardless of the two-dimensional radii used to extract 
the signal photons (e.g. Fukazawa et al. 1994; Paper I). In other words, the 2P modeling is 
likely to be more than a mere convention. 

Generally, the hot-phase and cool-phase temperatures, derived under the 2P assumption, 
are respectively close to T c and Th of eq.([TD obtain ed through the IP analysis. Furthermore, 
eq.flT]) holds for the 2P results as well (jlkebdl2002l ). The IP and 2P approaches are thus con- 
sistent in the sense that both yield essentially the same two characteristic temperatures, T c 
and Th. Nevertheless, their physical implications are much different. The IP situation would 
require a fine tuning between th e radiative cooling and th e putative heating mechanism, or 
"quasi-hydrostatic" gas cooling fjMasai fc Kitayamall2004l ). To realize the 2P condition, in 
contrast, the cool phase must be thermally insulat ed from the hot phase, and h eated in a sta- 
ble manner, so that it should neither evaporate (ITakahara fc Takaharalll979l ). nor collapse 
due to cooling. Yet ano ther possibility is th at the ICM temperature takes a range of values 
even at a single radius (IKaastra et al.ll2004l ); this is to be called "multi-phase" picture. 



In order to identify the CF-suppression mechanism, it is of basic importance to clarify 
which of the IP and 2P pictures (or the multiphase view) is closer to the reality. In spite of 
its importance, there have been few attempts to address this issue. We hence analyze in the 
present paper the XMM-Newton data of the Centaurus cluster. This object is best suited 
to our purpose, because it harbors the very prominent cool component and has extensively 
been studied with ASCA. Our strategy is to examine whether or not deprojected spectra of 
this cluster require the presence of multiple temperatures at each 3-dimensional radius. The 
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data will also allow us to reconfirm the lack of emission with temperature below the value 
of eq.©. 

In §2, we briefly summarize previous X-ray results on the Centaurus cluster, and describe 
the XMM-Newton observations in §3. Section 4, which forms the core of the present paper, 
is devoted to the description of data analysis and results. The results are discussed in §5, 
followed by a summary in §6. Throughout the paper, the Hubble constant is expressed as 
H Q = 72 km ~ l Mpc" 1 . The three-dimensional radius and the projected radius are denoted 
as R and r, respectively. Errors represent 90 % confidence limits unless stated otherwise. 



2. PREVIOUS X-RAY RESULTS ON THE CENTAURUS CLUSTER 



Located at a redshift of z = 0.0104 (ILucey. Currie fc Dickens Ill986l ; ISmith et al.ll2000i ) 
where 1' corresponds to 12.5 kpc, the Centaurus cluster (Abell 3526) is one of the most 
well studied objects among nearby clusters. It exhibits a roughly round shape, and its 
cD galaxy, NGC 4696, does not show significant nuclear activity. In X-rays, the Centauru s 

cluster has been observed exte nsively, wi th Uhuru (IGiacconi et al.lll972l ; iForman et al ll978h. 

OS08 dMushotzkv et al.lll978h. EE API jMitchell fc Mushotzk4l98(k Einstein jMatilskv. Jones fc Formal 



19851 ). EXOSAT ( lEdge fc Stewart! ll991a! Jbl). and Ginga ( lYamashitalll992l ). These observa- 
tions yielded an ICM temperature of ~ 4 keV, except in the central region (r < 6') where a 
very prominent cool component and an excess brightness are observed. 



Allen fc Fabian! (119941 ) studied the Centaurus cluster using the data from the ROSAT 



PSPC. They found some evidence of metallicity increase and excess absorption at the center. 
Through a deprojection analysis, they also derived 3-dimensional properties of the ICM under 
the IP modeling, reporting that the ICM is cooled down to ~ 1 keV toward the center. 



Using ASC A 



with the much improved energy resolution and the wider energy band, 



Fukazawa et al.l (119941 ) successfully resolved both the Fe-K line and low-energy emission 
lines, and detected hot and cool components with a temperature of ~ 4 and ~ 1 keV, 
respectively. The strong metallicity increase toward the center was reconfirmed. They also 
discovered that the 3-dimensional core region (R < 5') of the cluster is filled mostly with 
the hot component, whereas the cool component occupi es a m i nor f raction of the volume. 
Through a more detailed analysis of the ASCA data. Ilkebd (119951) showed that the IP 
density and temperature profiles derived with ROSAT (lAllen fc Fabian! 119941 ) underpredict 
the ha rd X-r a y flux from the central region, and hence fail to explain the 0.5-10 keV ASCA 
data. Ilkebd (119951 ) instead found that the ASCA spectra from various regions near the 



center can be described adequately by the 2P formalism employing the two components 
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identified by iFukazawa et al.l (119941 ) . This reinforced the presence of the hot phase in the 
three-dimensional core region, and justified the 2P approach. 

The 2P form alism on the Centaurus cluster was further refined and reinforced by 
Ikebe et al.l (119991 ). Incorporating two temperatures of ~ 4 keV and T c ~ 1.4 keV, 
together with a central metallicity increase, they successfully constructed a 3-dimensional 
ICM model that can simultaneously explain the R OSA T surface brightness and the ASCA 
annular spectra, namely the best- quality data sets ay ailable at that time. The excess ab- 
sorption at the center, clamed by iFabian et al.l (119941 ) based on the ASCA data, was not 
reconfirmed by llkebe et al.l (119991 ). 



Utilizing the unprecedented angular resolution of Chandra. ISanders &: Fabian! (120021 ) de- 
tected a "plume" like X-ray feature near the center of the Cen taurus cluster, reconfirming a 
previ ous observation with the ROSATHigh Resolution Imager (jSparks. Jedrzejewski Sz Macchetto 
1994J). Within 10" of the center including the "plu me" region, a cool plasma component with 
a temperature of ~ 0.7 keV was clearly detected (ISanders fc Fabianll2002l ). Employing the 
IP formalism, ISanders fc Fabian! (120021 ) also derived the ICM temperature distribution near 
the center (~ 3'), but the field of view of Chandra was not wide enough to determine the 
global temperature structure over the entire "cool" region (r < 6'). 

There are some indications that a subcluster centered on the elliptical galaxy NGC 4709 
is merging with the main body. Th is region exhibits an excess X-ray surface brightness and a 
temp erature increase to ~ 5 keV (jChurazov et al.lll999l ; iFurusho et al.l 1200 It iMolendi et al. 
20021 ). However, these effects, seen ~ 15' off the cen ter, are considere d to be small in regions 



which we analyze (r < 12'). In fact, using Suzaku, lOta et al.l (120071 ) showed that the ICM 
within the central ~ 12' can be regarded as hydrostatic: specifically, they set an upper limit 
of ~ 1,400 km s" 1 on any large-scale (~ 100 kpc) bulk motion in the IC M over this region, 
and r uled out a previously claimed detection of significant bulk motions (IDupke fc Bregman 
2006h . 



3. OBSERVATIONS AND DATA REDUCTION 



3.1. Observation 



The Centaurus cluster was observed with XMM-Newton on 2002 January 3 for a gross 



exposure of 45 ksec. The boresight was set to (a 



2000 £2000 



l2 h 48 m 49 s .3,-41°18'40".0). 



ters 


Striider et al. 


2001 


2001 


)• 
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The EPIC dataset obtained in the same observation were already analyzed by lMatsushita et al 



f )2007al ) mainly for metallicity distr i butio ns. The present paper emphasizes the temperature 
structure. Similarly, I Sanders et al.l (120081 ) utilized the RGS data from the same observation, 
as well as those from an additional 110 ks observation. 



3.2. EPIC Data Reduction 



3.2.1. Data processing 



The EPIC data were extracted using the standard SAS software of version 5.4.1. We 
selected events with pattern < 4 for PN and pattern < 12 for MOS, with which most of the 
valid X-ray events are accepted. According to instructions by the detector team, we further 
discarded events with bad flags, e.g., those out of the field of view and those detected at bad 
CCD pixels. 

Soft protons often cause the EPIC background to increase suddenly by up to two orders 
of magnitude. These sporadic "flares" must be removed carefully when analyzing extended 
sources. We accordingly produced 2-7 keV band light curves of the present data, excluding 
a central regio n (r < 8') where t he clu ster emission dominates. This particular energy band 
was shown by iKatayama et al.l (120041 ) to serve as a good measure of the EPIC background, 
and is roughly optimized in the present case to provide a high signal-to-ratio for proton 
flares; below ~ 2 keV the cluster signals become significant, while above ~ 7 keV the 
fluorescent backgro und lines in the PN spec tra b ecome a nu i sance . The results are shown in 
Figure [TJ Following IKatayama et al.l (120041 ) and iTakahashil (120041 ). we then discarded those 
time periods when the count rates deviate by more than 2 a from those in quiescent periods. 
After this screening, the net exposure time became 31 ksec with PN and 43 ksec with MOS. 

Several faint point sources are found in the field of view. We excluded regions where 
their signals are significant, typically within 15" of each source. 



3.2.2. Background 

Since the cluster emission fills nearly the entire EPIC field of view, background spectra 
must be constructed using some other "blank sky" data. For this purpose, we utilized the 
observation of PG 1115+080 for PN (Obsld=0082340101), and that of Vik 59 for MOS (Ob- 
sld=0107860501), because these data sets showed a close similarity to the present Centaurus 
data in terms of high-energy (11-14 keV for PN and 10-12 keV for MOS) count rates accu- 



-7- 



mulated over peripheral regions of the fields of view. 



After removing proton flares and celestial point X-ray sources in the same manner, we 
extracted background events from these data sets, using the same regions on the detectors 
as employed to accumulate the on-source data. The background derived in this way is 
estimat ed to be accurate gen erally to within ~ 3%, except at some particular spectral 
regions (IKatayama et al.ll2004l ). as detailed in §3.2.41 where we employ 8% . This accuracy is 
sufficient for our purpose, because the cluster emission is much brighter than the background, 
e.g., by more than a factor of five, at least up to ~ 5' from the aimpoint in energies below 5 
keV. Thus, the background uncertainty therein is at most a few percent of the cluster signal. 



3.2.3. Deprojection 



To discriminate the IP and 2P conditions, it is essential to remove the projection effects, 
using so-called deprojection procedure . To create deprojected spectra, w e employed the 
standard "onion peeling" method (e.g. Ilkebe. Bohringer fc Kitayamal 120041 ). assuming that 
the object is spherically symmetric, and the plasma is uniformly distributed in each thin 
shell region. No model spectra were assumed in our deprojection procedure. Specifically, at 
each energy E, we calculated the deprojected spectrum Sj(E) in the j-th shell as 



A? 



S J (E) = J2v j , k A k (E) 



(2) 



where Ak(E) is the projected spectrum accumulated over the k-th annulus, {T>j^} is a 
triangular matrix derived by inverting the projection matrix, and N denotes the outermost 
annulus. 

Below, we utilize the data within 12' (150 kpc) of the X-ray centroid at (12 h 48 m 49 s .2, -41°18'44") 
which coincides with the NGC 4696 nucleus. The emission outside 12', Sn(E), was calcu- 
lated assuming that the plasma is isothermal and isoabundance at R > 10', and that the 
ICM density profile is expressed by a j3 mode l determined joint ly with ASCA and ROSAT, 



namely j3 = 0.57 and the core radius of 7'. 3 (Ilkebe et al. 
produce the XMM-Newton images themselves (jTakahashi 
{T>ij} are numerically given in Appendix A. 



1999} . These paramet 



ers can re- 



20041 1 . Actual matrix elements of 
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3.2.4- Systematic Errors 



In order to take int o account calibration uncertainties, conservative systematic errors 
of 2% fjLumb et al.ll2002l ) were assigned to the on-source spectra in energies above 0.6 keV. 
Below that energy, the systematic error was increased to 5% t o represent uncertaint ies in 



the quantum efficiency around the oxygen absorption edge (e.g. iKatayama et al.l 120041 ) 



To further consider systematic uncertainties involved in the background subtraction 
process, we assigned to the background data the systematic errors as given in Table HJ The 
errors are thus taken to be 8% ( §3.2.21) in those energy bands where the background is 
dominated by instrumental fluorescence lines, because their intensities are known to vary 
rather independently of that of the background continuum. In energies below ~ 1 keV, the 
background spectra are significantly contributed by diffuse soft X-ray background, of which 
the brightness varies from sky to sky. There, we hence employed the same 8% systematic 
error, although the background in this energy range is generally negligible compared to the 
bright cluster emission. In other energy regions, we employed a nominal value of 3% ( §3.2.2p . 



3.3. RGS Data Reduction 

Again using the SAS software, we processed the RGS data, and extracted the first and 
second order spectra within r = 60" (in the cross-dispersion direction) of the center. The 
spectral range of 6-23 A was analyzed, in order to utilize many atomic lines therein. 

The RGS backg round was derived from blank sky fields which were made available by 



the instrument team (ITamura. den Herder fc Gonzares-Riestrall2003l ). Because these authors 



reported that the rms fluctuation of the RGS background is not more than 30% over the full 
spectral range, we assigned a systematic error of 30% to the RGS background. In reality, 
the emission from the Centaurus cluster exceeds the background by more than an order of 
magnitude over the interested spectral range; therefore, the background uncertainties have 
little effects on the results. 
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4. DATA ANALYSIS AND RESULTS 
4.1. EPIC Data Analysis 

4.1.1. EPIC Spectra 

Figure [2] shows the PN spectra of the Centaurus cluster, obtained by accumulating the 
on-source events over a series of concentric annular regions centered on the NGC 4696 nucleus 
( §3.2. ip . and then subtracting the background described in §3.2.21 The radii defining these 
annuli are given in the figure caption. The figure also shows their ratios to those averaged 
over an outer region of 5' < r < 12', where the cool component is hardly seen. The spectra 
from these regions do not bear strong atomic lines except the Fe-K line at ~ 6.7 keV, 
and their spectral shapes are very similar to one another. Therefore, these regions at r > 5' 
must be filled with an approximately isothermal hot plasma. The figure also confirms correct 
background subtraction, because wrong subtraction would make the hardest-band (> 7 keV) 
spectra significantly depend on the radius. 

Toward the center, in contrast, we can see prominent increases in the line intensities; 
the Fe-K line enhancement is caused by the metallicity increase, while those of Si-K (at 
rest-frame energies of 1.86 and 2.01 keV), S-K (2.46 and 2.62 keV), and Fe-L (~ 1 keV) lines 
are mainly due to the emergence of the cool component. 

At the very central region (r < 1'), hereafter called "cD region", Fe-L lines with energies 
below 1 keV become prominent, indicating the presence of a still cooler plasma with a 
temperature of < 1 keV. The same inference can be obtained from the drastic increases in 
the H-like to He-like intensity ratios of the Si-K and S-K lines. We therefore qualitatively 
reconfirm t he Chandra detection o f the coolest (~ 0.7 keV in temperature) component at 



the center ( iSanders fc Fabianl 120021 ) 



Applying the procedure described in §3.2.31 to the spectra in Figure (2], we derived de- 
projected spectra from 11 concentric shells. Some of them are presented in Figure [21 in 
comparison with the non-deprojected spectra. The deprojection obviously degrades the data 
statistics, but does not alter the spectral shape very much; the only eye-catching difference 
is the reduction in the hard X-ray (> 3 keV) flux of the innermost spectra (denoted [cd]), 
caused by the removal of contributions from the foreground/background hot emission. 
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4-1-2. Single-phase analysis 



As first-cut evaluation of the annular spectra of Figure HJ we conducted a conventional 
single-temperature analysis. This is meant to approximate IP conditions, but only crudely, 
because three effects, coupled with the temperature and abundance gradients, would make 
each spectrum deviate from an exact isothermality; the most dominant one is the projection 
effects, while the less important ones are the finite angular resolution of the X-ray telescope 
which mixes up signals from adjacent shells, and the finite thickness of each annulus. Thus, 
the single-temperature representation should be regarded as a crude approximation here. 

We fitted the three projected EPIC spectra from each annulus with a single-temperature 
APEC model fjSmith et al.ll200ll ) with photoelectric absorption, simultaneously but allowing 
the model normalization to differ among the three detectors. The abundances of O, Ne, Mg, 
Si, S, Ar, Ca, Fe, and Ni wer e all le ft free, while those of C and N were fixed at the solar 
values by lAnders fe Grevessd (119891 ). The redshift was fixed at 0.0104, and the hydrogen 
column density at the Galactic-line-ofsight value, 8.1 x 10 20 cm -2 , because iMatsushit a et al. 



(j2007al ) found the absorption to be grossly consistent with it. Since some systematic residuals 
were observed around the Fe-K line region, we incorporated small gain correction factors, 
+0.3% to PN, +0.6% to MOS1, and +0.3% to MOS2, all within the reported calibration 
uncertainties. These improved the fit chi-squared, without significantly affecting the best-fit 
parameters. 

As expected from Figure fits to the cD region spectra left large residuals around the 
Fe-L energies. We therefore added, in these cases, another cool APEC component with a 
temperature (left free) of ~ 0.7 keV, of which the abundances are tied to that of the major 
component. As a result , the fits were much improved, in agreement with the report by 
Sanders fe Fabian! (120021 ). This modeling should be distinguished from the 2P approach to 
be examined later. 

The IP model (plus the 0.7 keV component) reproduced the spectra moderately well, 
although the fits are not necessarily acceptable. Figure H] (filed circles) shows radial profiles 
of the obtained parameters. Thus, the plasma temperature approaches a constant value 
of ~ 3.8 keV toward outer regions, while drops i 
from Figure [2b and with the previous results (e.g 



Fabian et al. 


1994; 


Fukazawa et al. 


1994 



Ikebe et al.lll999t iFurusho et al.ll200ll ). As mentioned above, the ~ 0.7 keV plasma compo- 



nent was needed at r < 1'. 

In Figure HJ the abun dances of Fe and Si are seen to increase toward the center as 
first revealed with ASCA ( jFukazawa et al.lll994l ). while that of oxygen is spatially rather 
constant. These results are consistent with the single-temperature analysis results from 
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Matsushita et al.l (l2007al ). Although ISanders fc Fabian! (120021 ) reported a Chandra detection 



of a metallicity decrease within r ~ 45" (down to ~ 0.4 solar at r < 10"), such an effect, if 
any, is less significant in the EPIC data. 

We performed the same model fitting to the deprojected spectra. Since the spectra are 
now free from the projection effects, the single-temperature modeling is equivalent to IP 
conditions, assuming that each deprojected shell is thin enough and the signal mixing due 
to the telescope response is negligible. This analysis yielded the data points represented in 
Figure H] by open diamonds. Thus, the overall results are very similar to those obtained 
with the annular spectra. The fits became generally better, due to the reduced statistics 
of the deprojected spectra, and to the removal of fore- and background contributions. The 
latter also makes the deprojected temperatures systematically lower (typically by 0.3-0.5 
keV) than the non-deprojected values. 



4-1.3. Two-phase analysis 



As an alternative approach, we fitted the EPIC spectra by a 2P model, i.e., a sum of 
two APEC components with different temperatures, plus the ~ 0.7 keV component at < V . 
For each annulus, the cool component temperature T c was left free, while the abundance of 
each element was assumed to be the same among the components. Like in the IP analy- 
sis ( §4.1.2p . the C and N abundances were fixed; so were the redshift and the photoelectric 
absorption. When we first let the hot component temperature Th float as well, it became 
rather unconstrained in some out er regions, and scattered over a range of ~ 3 to ~ 7 keV. 
This is due to coupling with T c ([Matsushita et al.ll2007al ). However, the value was gener- 
ally consistent, within rather large errors, with 3.8-4.0 keV obtained in past observations. 
We hence fixed it at Th = 3.8 keV after the ASCA results (jlkebe et al.lll999l ) and our IP 
fits (Figure U). 

By this 2P model, both the non-deprojected and deprojected spectra were reproduced 
reasonably well in all regions. The obtained parameters, shown in Fig ure are generally 
consistent with the two-temperature results by iMatsushita et al.l (|2007aJ) within some differ- 
ences in the modeling. Thus, the fit goodness is comparable to, or slightly better than, those 
obtained with the IP analysis. The derived abundance profiles are essentially the same as the 
IP results, and th e radial profile of Fe (Fig. [5b) agrees very well with the ASCA 2P results 
(jlkebe et al.lll999l ). The temperature differences between the projected and non-deprojected 
spectra became smaller than in the IP analysis, and often insignificant within errors. This 
is consistent with a 2P condition, in which an annulus and the corresponding deprojected 
shell are expected to share nearly the same two temperatures. 
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Figure[6ti shows radial profiles of the emission measures (per unit volume) of the hot and 
cool components, denoted Qh and Q c , respectively. Thus, the cool component is required 
within 6', where the IP temperature decreases significantly in our IP result (Figure H]). 
The derived T c stays at ~ 2 keV over 2' — 6' (though rather poorly determined outside 3'; 
Figure [5^), while it decreases down to T c ~ 1.6 keV at the center. The la tter agrees with 
the c entral IP temperature from our IP fits (Figure H]) and from ROSAT (lAllen fc Fabian 
19941 ) . At the cD region, the third 0.7 keV component is again required. In contrast, neither 
the ring-sorted spectra nor the deprojected thin-shell ones required the cool component in 
regions outside 6'. Accordingly, Figure [5] outside 6' simply reproduces the IP results. 

I ncidentally, the 2P an alyses of the ASCA data yielded some what lower values of T c ; ~ 1 
keV (IFukazawa et al.lll994l ). or 1.4 ± 0.2 keV (llkebe et al.lll999l ). This is partially because 
ASCA was not able to separately detect the 0.7 keV component. The remaining difference 
may be due to the different plasma codes ([Matsushita et al.ll2007a ). the MEKAL and APEC 
codes used in the ASCA study and the present analysis, respectively. 



To characterize the obtained 2P solution, let us introduce, after llkebe et al.l (119991 ). the 
volume filling factor r] c of the cool phase. Using the emission measures, this quantity is 
defined as 

Qc = £n 2 c r] c , Q h = £nl(l - r) c ) , (3) 

where and n c are the densities of the hot and cool phases, respectively, while £ is a 
factor o f order unity reflecting the metallicity. Assuming a pressure balance between the two 
phases (IFukazawa et al.lll994l ). namely, 



n c T c 



n h T h 



(4) 



we can calculate rj c as 



i + i£V7£ 



Qc 



(5) 



Figure [6)3 shows the radial profile of rj c , thus calculated using eq.©; in deriving this 
result, we fixed T r at 2.0 keV in shell regions outs ide R = V . The figure reconfirms the ASCA 
results ( IFukazawa et al .11 19941 ; llkebe et al.lll999l ). that the hot phase dominates the volume 
down to a radius rather close (R ~ 1'.5 or 18 kpc) to the center. The filling factor implied 
by the XMM-Newton data is higher by a factor of ~ 3 than the ASCA determination; as can 
be understood via eq. (JSJ), this is due partially to the higher value of T c . 
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4.2. Analysis of Thick-Shell Spectra 

As we have seen in §4.11 the two descriptions, IP and 2P, are consistent with each other 
in the sense that both find essentially the same set of three plasma temperatures, ~ 3.8 keV, 
1.7-2.0 keV, and ~ 0.7 keV. Outside ~ 6' (~ 75 kpc), the ICM can be approximated as 
isothermal at ~ 3.8 keV. Within ~ 6' toward the cluster core, the ICM becomes gradually 
cooler down to ~ 2 keV (if employing the IP description), or the ~ 2 keV plasma becomes 
intermixed with the hot 3.8 keV component with a progressively larger r) c (if employing the 
2P description). In the cD region, both scenarios require the additional 0.7 keV plasma. 
Up to this stage, we are unable to tell which of the two approaches is more favored by the 
data. Although the 2P fit to the deprojected spectra often gives a slightly lower reduced chi- 
squared than the IP fit (FigureHb vs. Figure[5b), the difference is not necessarily statistically 
significant, and could be due to the finite temperature gradient within each shell. 



4-2.1. Construction of thick- shell spectra 

In order to better distinguish the IP and 2P modelings under higher statistics, we took 
the deprojected spectra from consecutive five thin shells covering altogether R — V — 5' ([e] 
though [i] in Figure [2]), which have IP temperatures of 2.03, 2.48, 2.65, 3.10, and 3.30 keV 
respectively ( §4. 1 .2[) . We then summed them into a single spectrum which now represents a 
thicker 3-dimensional shell in the cool core region; this is hereafter called "Shell C" . The cD 
region, together with the "plume" structure, is excluded to avoid complexity introduced by 
the 0.7 keV component. Since the statistical errors associated with the adjacent thin-shell 
spectra are not mutually independent, we estimated the errors to be a ssigned to the th ick- 



shell spectra by properly considering error propagation utilizing eq. (121) (ITakahashil 120041 ) . 



The Shell C spectra, one for each detector, should contain the 3.8 keV component if 
the 2P scenario is correct. If instead the IP picture is more appropriate, the Shell C spectra 
would rule out contributions from such a high temperature component, because the IP 
temperature changes across Shell C from 2.0 keV to 3.3 keV but not higher (Figure H]). We 
expect that the improved data statistics, achieved by the spectral summation, allow us to 
distinguish these two cases. 

In addition to the improved statistics, analyzing a thick shell has two more advantages: 
the signal mixing caused by the finite telescope resolution is reduced, and the errors due to 
spectral changes within individual thin shells, if any, nearly cancel out. As an instructive 
exercise of the latter effect, let us assume, for example, the 3rd and 4th annuli constituting 
Shell C to have a difference by 5W in their Fe-K line equivalent widths, and examine what 
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would happen if we artificially assigned the two annuli with a same value W that is the 
emission-measure-weighted average between the two shells. This exaggerates the procedure 
of averaging out a radial spectral change over a finite shell thickness. The errors caused 
by the averaging are expected to propagate through deprojection into the 1st through 4th 
thin shells. However, as shown in Appendix A, a mathematical estimate using actual matrix 
elements of eq.([2]) confirms that the errors in the four thin shells roughly cancel out, and 
hence the final thick Shell-C spectrum suffers in the present case no larger errors in the line 
equivalent width than ~ 0.015W 7 . Since the same estimate applies to a temperature gradient, 
which is at most ~ 25% in each thin shell (Figure H|) over 1' — 5', the residual temperature 
errors in the Shell C spectrum are estimated to be at most ~ 0.3%, and hence negligible. 



4-2.2. Single-phase modeling 

Let us examine the Shell C spectra against the IP view. We hence took the best-fit 
single-temperature APEC model for each of the 5 constituent thin shells, and convolved it 
with a response of each detector that is weighted by the surface brightness of the projection 
of the relevant thin shell. In the same way as the actual thick-shell data, we then summed up 
the 5 simulated thin-shell spectra of each detector into a single one, to be called "synthetic 
IP spectrum". 

Figure [7^ directly compares these synthetic IP spectra with the actual Shell C spectra of 
the corresponding detectors. Thus, we observe a moderately good agreement, but the actual 
Shell C spectra exhibit some excess above the synthetic IP model toward higher energies. In 
fact, their difference, in terms of chi-squared summed over the three detectors (PN, MOS1 
and MOS2), becomes x 2 / u = 739/378 as summarized in Table [2J 

The data vs. model comparison in Figure [T^i incorporates no adjustment, since the syn- 
thetic IP spectra are uniquely specified by the best-fit models describing the 5 constituent 
thin shells ( §4.1.21) . However, the agreement might be improved by adjusting the model 
parameters. We therefore allowed to vary freely the normalizations of the five APEC com- 
ponents which constitute the synthetic IP spectrum. Then, the agreement (now to be called 
"a fit") was improved to x 2 / u = 724/373 (Table [2]); the normalizations of the five compo- 
nents changed by a factor of 1.65, 0.25, 0.84, 1.10, and 1.26 in the increasing order of the 
temperature, while the summed normalization remained unchanged within ~ 2%. Thus, the 
coolest and hottest of the five components tend to be enhanced, implying a 2P-like condition. 
To retain a stable fit convergence, we did not attempt to let the 5 temperatures float. 
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4-2.3. Two-phase modeling 

In order to next examine the 2P approach, we fitted the actual Shell C spectra with a 
sum of two APEC models, exactly in the same manner as in §4.1.3t the two temperatures 
and the two normalizations were left free, as were the abundances of O, Ne, Mg, Si, S, Ar, 
Ca, Fe, and Ni which are assumed to be common to the two components. The relative 
normalizations among the 3 detectors were left free as before. As a result, the degree of 
freedom decreased by 15 compared to the synthetic IP modeling. The two components were 
convolved with the same response, weighted by the X-ray surface brightness. 

As shown in Figure [7b and Table [2J the obtained fit (x 2 /V = 688/363) is significantly 
better than that with the synthetic IP model, even considering the reduced degree of freedom. 
In fact, the high energy excess, observed in Figure [TJi, has been reduced in Figure [7b. The 
two temperatures have been obtained as Th = 4.02^33 keV and T c = 2.06^q;?q keV. While 
these values agree well with the 2P results from the individual thin shells (Figure E]), the 
obtained Th is significantly higher than the temperature range (up to 3.30 keV, with a typical 
error of ±0.15 keV) involved in the IP description of Shell C (Figured]). Although the error 
ranges associated with the derived Th and T c are somewhat (typically by ~ 40%) under- 
estimated due to rather large values of the reduced chi-squared (~ 1.9), the above inference 
remains valid even if the error ranges are enlarged by ~ 40%. 

Just for consistency, we applied in the same way the 2P fit to the synthetic IP spectra, 
to which Poissonian noise was added. This yielded T c = 2.3 keV and Th = 3.3 keV, although 
the associated errors are rather difficult to estimate. These two temperatures are both 
contained within the IP temperature distribution across Shell C (Figured]), from 2.03 keV 
to 3.30 keV, that we employed in constructing the synthetic IP spectra. This difference 
between the actual Shell C spectrum and the sythetic IP spectrum suggests that the hot 
3.8 keV component was already present in the 5 thin-shell spectra, and has become more 
significant by the spectral summation. 

In the above 2P fit to the actual Shell C spectra, the two components were treated to 
have the same abundances, and were convolved with the same response. However, the true 2P 
condition would be somewhat different, since the cool component, more weighted toward the 
metal-enriched inner regions than the hot one, must have higher abundances when averaged 
over Shell C. To better express this condition, we returned to the original best-fit 2P models 
to the five thin-shell spectra, convolved them with projected responses of the respective thin 
shells, and summed the results over the 5 shells. These simulated data, involving two fixed 
temperatures (2.0 and 3.8 keV) and properly considering the radial abundance changes, are 
to be called "synthetic 2P spectra", because they are constructed exactly in the same way 
as their IP counterparts. Although these synthetic 2P spectra do not involve any adjustable 
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free parameter, they have reproduce the Shell C spectra with x 2 jv = 678/378, even better 
than the above 2P fit, and also than the synthetic IP modeling. 

4-2.4- Error renormaliztion 

From the statistical viewpoint, even the synthetic 2P simulation obtained in § 14.2.31 
does not give an acceptable reproduction of the data (Table [2]). This hampers correct error 
evaluations, and makes it difficult to quantitatively compare different modelings. In Fig- 
ure U\ fit residuals are observed around strong emission lines, such as Si-K, S-K, and Ar-K. 
Therefore, the inadequate fit goodness is likely to be caused mainly by subtle calibration 
inaccuracy in reproducing the detailed line profiles, rather than model inappropriateness. 
Under the higher data statistics, the calibration uncertainties presumably exceeded the sys- 
tematic errors of 2% already assigned in §3.2.41 Although our data analysis utilized rather 
old SAS version (5.4.1), a reanalysis of a few representative spectra from the central region, 
using an updated version 7.0.0 SAS software (processing and responses), showed rather in- 
significant improvements of these problems. Furthermore, the fitting results derived with 
the two versions agreed within errors. Accordingly, we retain the version 5.4.1 analysis, and 
assigned 5% of the source counts as a new systematic errors to the whole energy range, and 
repeated the overall analyses. This procedure is meant to make the fit acceptable when us- 
ing the model that is best preferred by the data, and to make securer the distinction among 
different models. 

By this error renormalization, the values of x 2 have reduced appreciably, while the radial 
profiles of various plasma parameters remained unchanged. Specifically, the comparison of 
the Shell C spectra with the remade synthetic spectra gives x 2 = 418 and 381 for IP and 2P, 
respectively, both with v = 378 (Table [2]). The latter is statistically acceptable, and gives a 
chi-squared which is smaller by 30 than the former. Thus, the 2P modeling of the Shell C 
spectra is statistically superior to the IP one. 

4-2.5. Thick Shell P spectra 

In a similar way, we derived EPIC spectra from another thick shell, called Shell P, which 
sums up three consecutive thin shells covering R — 3' — 6'. This Shell P is peripheral to the 
central cool region. The summed Shell P spectra were then compared with the corresponding 
synthetic IP and 2P spectra, constructed in the same manner as before from the best-fit 
solutions to the 3 constituent thin shells. The comparison of the actual spectra with the 
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synthetic IP and 2P spectra gave x 2 = 511 and 496 (without the error renormalization), 
respectively, both with v = 378. The difference in chi-squared, though smaller than in the 
Shell C case, amounts to 15, and remains ~ 11 even when the systematic error is renormalized 
as in 

We also fitted the Shell P spectra with the 2P model, in the same way as for the Shell C 
case but fixing at 3.8 keV as we did in t HT51 We then obtained T c = 2.10^ keV, 
together with x 2 / u — 494/364. This value of T c again agrees with those obtained by the 
2P fits to the constituent thin-shell spectra (Figure [5]), while it is significantly below the 
distribution of the IP temperature across Shell P, 3.0-3.5 keV (Figure H]). Therefore, the 
cool (~ 2.0 keV) component is inferred to be present even in this relatively outer region. As 
a consistency check, we fitted the synthetic IP spectra in Shell P with the same 2P model, 
and obtained T c = 3.2 keV which agrees with the input IP temperature range. 



4-2.6. Possible artifacts 



We have so far found that the two sets of thick-shell spectra both prefer the 2P view to 
the IP scenario. However, the results could be subject to various artifacts. 

An immediate suspect is that the appearance of the hot component in Shell C is an 
artifact caused by insufficient background subtraction. However, this would not explain the 
presence of the cool component in Shell P. Furthermore, the constant spectral shape in r > 5' 
( §4.1.1( Figure [2]) argues against this possibility. For a more quantitative examination, we 
remade the deprojected spectra with the background artificially under-subtracted, by 3% 
which is comparable to the background systematic errors estimated in §3.21 In fact, this 
little affected the results obtained so far: a comparison of the remade Shell C spectra with 
the 2P synthetic spectra ( §4.2.31) gives nearly the same x 2 as before, changing only by ~ 2. 
Thus, our results are robust against the background uncertainty. 

Another concern is possible deviations from the so-far assumed spherically symme- 
try. Indeed, the ICM te mperature of the Centaurus cluster exh i bits significant nqn-axi 



symmetric distribu tions (IFurusho et al.1 1200 ll ; iMolendi et all |2002| ; ISanders fc Fabian! 12002 



Fabian et al.ll2005l ). Then, even if the ICM were in a IT condition, Shell C would sample a 
wider range of temperatures, to mimic a 2T condition. However, as far as we consider an an- 
nular region of 4' < r < 8' that is crucial to this issue, th e ICM temperature (approximated 
as IT) is azimuthally constant typically within ±15% (IFurusho et al.ll200ll ). To examine 
this issue using the present data, we derived a temperature map, presented in Figure 
by analyzing the projected EPIC spectra in four sectors by a single-temperature model. 
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It reconfirms a mild temperature decrease from northwest to southeast directions, with an 
amplitude comparable to those measured previously. According to a simple estimation, a 
temperature variation of this amount is insufficient, by at least a factor of two, to explain 
the difference between the two values of T h = 4.02^33 keV and T h = 3.03 keV, obtained 
via the 2T fit to the Shell C spectrum and its synthetic IT data, respectively. Thus, the 
preference for the 2T view remains intact. 

The finite angular resolution of the XMM-Newton X-ray telescope could also blur the 
Shell C boundary, and cause contaminations by signals from the outer hotter region. How- 
ever, Shell C is much thicker (4') than the blurring width (~ 15"), and the surface brightness 
decreases outwards. As a result, the contamination of the outer hotter emission into Shell C 
is estimated to be at most ~ 2%, which is within the renormalized systematic uncertainty 
of 5% introduced in § 14.2.41 

In the deprojection process (§ I3.2.3p . we assumed the region outside 12' (150 kpc) to 
have a constant spectrum that is identical to that measured just inside it. This is based 
on the fact that the azimuthally-averaged ICM temperature of this cluster, measured with 
the ASCA GIS utilizing its wide field of view and lo w background, is essentially constant 



at ~ 3.9 keV from r = 14' up to r = 40' (0.5 Mpc; llkebe et al.lll999l ). Of course, devia- 
tions from this assumption would in principle affect, via eq.(J2]), the deprojected spectrum 
in every shell. However, our numerical estimate in Appendix A indicates that the emission 
from the R > 12' region contributes only ~ 1% to the Shell C spectrum, and hence any 
realistic error associated to it is considered negligible. Furthermore, the ICM temperature 
has been confirmed in many clusters to decrease beyond ~ 0.2 t imes the virial radius (e,g 



Markevitchlll996l : iDe Grandi fc Molendil 12002 ; ISato et all 12007b : iLeccardi &: Molendil 120081 ) 
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or r > 0.5 Mpc in the present case. Then, we are likely to be overestimating the outermost 
ICM temperature, and hence over-subtracting the hot (3.8 keV) emission as foreground and 
background contributions. This makes our argument simply more conservative. 

A more difficult issue is asymmetry along the line of sight. If the object has a prolate 
shape in the depth direction, the deprojected shell C spectra would sample effectively outer 
(hence hotter) regions, and would involve the 3.8 keV component even if it is in a IP 
condition. Similarly, an oblate IP condition would make Shell P sample the cool emission. 
However, neither case can explain the two thick shells simultaneously. 

Finally in § 14.2.11 possible spectral changes within each thin shell was already shown 
to become negligible through the construction of thick shells. As an overall confirmation to 
address this issue, in Appendix B we have numerically simulated a IP cluster, and analyzed 
the fake data in the same manner as for the actual data. Then, the simulated Shell C spectra 
have been reproduced by the IP view, without demanding the 2P formalism. As a logical 
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contraposition, the preference of the 2P modeling by the actual data demands the actual 
object not to be in the IP condition. 

From these examinations, we conclude that the data favor the presence of the hot 
~ 3.8 keV component in the 3-dimensional core region (R = V — 5'), and the cool ~ 2 keV 
component in the peripheral region (R = 3' — 6'). In other word, the 2P picture is considered 
to describes the actual physical state better than the IP view. 



4.3. Excess Absorption 



At the center of the Centaurus cluster, some previous works reported the presence of 
excess X-ray absorpti on, up to a few times ~ 10 21 cm" 1 (depending; on the modeling;), above 
the Galactic column (lAllen fc Fabianlll994l ; iFabian et al.lll994l ; ISanders fc Fabianll2002l ); this 
was taken as evidence for the CF. However, the effect was not confirmed in other works 



( llkebe et al.lll999l ). Likewise, we have here experienced no problems in our spectral analysis 
with the absorption fixed at the Galactic value. 

Generally, a somewhat higher absorption is often derived when a tw o- (or multi-) tem- 
perature emission is fitted with a IP model (e.g.. lMatsushita et al.ll2007al ). However, some of 
the previously reported excess column densities in the Centaurus cluster considerably exceed 
such modeling uncertainties. 

We examined this issue using the non-deprojected and deprojected spectra from regions 
within 10" and 30" of the center, respectively. Specifically, we refitted these EPIC data using 
the 2P plus 0.7 keV model, applying an intrinsic absorption factor to both components in 
addition to that due to the Galactic column. (The hot component can be neglect here.) As 
a result, the fits were little improved in either case. The obtained excess column density 
in the rest frame is 3.2+ 2 ;° x 10 20 cm 2 for the non-deprojected spectra within 10", and 
0.51q 5 x 10 20 cm -2 for the deprojected spectra within 30". These are significantly lower than 
the previously claimed values, > 1 x 10 21 cm -2 , and the case of the deprojected spectrum is 
consistent with no excess absorption. Then, what is responsible for this discrepancy? 

Regardless of the modeling (IP or 2P), the XMM-Newton data reveal that oxygen in 
the central region is deficient relative to iron and silicon ( §4.1.2} §4.1.31) . Nevertheless, the 
previous works usually assumed metals in the ICM to obey the solar abundance ratios, simply 
because of inadequate energy resolution. This presumably caused data deficits around the 
oxygen line (~ 0.65 keV) relative to the model predictions, and forced previous investigators 
to introduce an artificial excess absorption in an attempt to suppress the over-predicted O-K 
line flux in their models. Actually, if we constrain the ICM to follow the solar abundance 
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ratios, the present data clearly require an excess absorption by ~ 1.4 x 10 21 cm -2 , although 
the fit becomes worse (with x 2 / v — 1-38) compared to the case of the Galactic absorption and 
free abundance ratios (x 2 /^ = 1-13). Thus, the reported excess absorption in the Centaurus 
cluster is considered an artifact, caused mainly by the non-solar abundance ratios around 
the cD galaxy. 



4.4. RGS Results 

Because the surface brightness of the Centaurus cluster is strongly peaked at the center, 
we can utilized the RGS data to better constrain the temperature structure of the ICM in the 
cD region, although the available spatial information is limited to 1-dimensional projections. 
We extracted the first and second order spectra from RGS1 and RGS2, over a strip of 2' in 
width centered on the NGC 4696 nucleus. The cross-dispersion direction is at a position angle 
of ~ 20°, so that the data accumulation strip partially covers the plume. The background 
was subtracted as described in §3.31 

Figure [9^ shows the obtained RGS spectrum, which combines the 1st and 2nd order 
spectra from the two RGS units. Although the energy resolution is somewhat degraded due 
to the finite angular extent of the source, we clearly observe many atomic emission lines. In 
particular, the spectrum bears a strong K-line at ~ 19 A from H-like oxygen ions (OVIII), 
but lacks those from He-like ones (OVII) which would be emitted strongly at about 22 A by 
plasmas with a temperature lower than 0.3 keV. This indicates that the cluster core region 
is devoid of such very cool plasmas, contrary to the prediction by the CF hypothesis. 

Since we already know that at least two plasma temperatures (0.7 keV and 1.7-2.0 
keV) are required to reproduce the emission from the cD region, we fitted the RGS spectra 
with a simple 2P plasma emission model; we hereafter call this particular modeling "quasi- 
2P fit", because it indeed invokes two APEC components but the employed temperatures 
are different from those composing the 2P model employed in §4.1.31 and §4.21 The four 
spectra (1st- and 2nd-order spectra from RGS1 and RGS2) were prepared separately, and 
fitted jointly. We neglected for the moment the hot (3.8 keV) component, because its high 
temperature and large angular extent make low-energy lines rather weak, and hence its 
parameters are difficult to constrain with the RGS data. The abundances of O, Ne, Mg, Si, 
Fe, and Ni were le ft free, but const rained to be common between the two components. After 



the EPIC results ( ITakahashil 120041 ) . we fixed the S, Ar, and Ca abundances at 1.7, 1.5, and 
2.4 solar, respectively, because K-lines of these elements are outside the RGS wavelength 
range. The absorption is assumed to be Galactic, as confirmed in §4.31 To deal with the 
spatial extent, we blurred each spectral component with a Gaussian of a free width a. 
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The quasi-2P fit to the RGS spectra is shown in panels (b) and (c) of Figure [9], and 
the obtained parameters are summarized in Table [31 Thus, the overall emission-line features 
have been reproduced to a reasonable extent. Although the fit is not yet fully acceptable, 
this is mainly due to some unessential factors, such as our simple spectral blurring, the 
projection effects, and calibration uncertainties. 

The obtained two temperatures, 0.8 keV and 1.7 keV, agree well with the lowest two 
characteristic values measured with the EPIC from the cD region. The Gaussian widths 
of the 0.8 keV and 1.7 keV components, ~ 20" and ~ 65" respectively, can be understood 
as reflecting the angular extent of the ~ 0.7 keV and ~ 1.7 keV component detected with 
the EPIC. The RGS-determined abundances of Mg, Si, and Fe are lower by a factor of ~ 2 
than the values obtained with the EPIC data. The discrepancy is at least partially relaxed 
by considering the hot (3.8 keV) component, which does not emit prominent lines, but is 
estimated to contribute to the continuum up to 30% and reduce equivalent widths of the 
lines by that amount. T he remaining abundance discrepancy could be due to the abundance 



drop at the very center (ISanders fe Fabianll2002l ). In any case, the temperatures, determined 



primarily by the line intensity ratios, are not affected. 

The RGS spectra may be examined for the presence of other temperature components, 
besides the ~ 0.8 keV and 1.7 keV ones constituting the quasi-2P fit. The RGS s pectr a 
were therefore analyzed with the multi-temperature fit procedure of 
whi ch is a kind of differe ntial emission measure analysis like those in 



Tamura et all (120031) 



Peterson et all fl2003h 



and iKaastra et al.l (120041 ) . Specifically, we prepared seven plasma emission components of 
which the temperatures are given as T , 1.5T , (1.5) 2 T , .., and (1.5) 6 T , with T being the 
lowest temperature. We then allowed to vary the normalizations of the seven components, 
as well as T to which the remaining 6 temperatures scale. The abundances of O, Ne, Mg, 
Si, Fe, and Ni were left free, but constrained to be common among the seven components. 
Each spectral component was blurred with a Gaussian, of which the width is left free but 
constrained to be less than 66" as obtained in the previous fit (Table [3]). 

As shown in Figure [101 the best fit model obtained from this analysis is physically 
very close to that obtained with the quasi-2P model. The fit has been improved little 
(x 2 /^ = 976/718; cf. Table [3]) by considering extra temperature components, with only 
the two components (~ 0.8 keV and 1.7 keV) remaining significant. Note that these two 
temperatures are constrained to have a ratio of 1.5 2 = 2.25, while their absolute values are 
left free. The 1.1 keV component has turned out to be weak in contrast to the two adjacent 
ones, implying that the significantly detected two components are discrete entities rather 
than representing a continuous temperature distribution. The 2.55 keV component is not 
significant, either. 
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The RGS results constrain the 3.8 (=1.7 x 1.5 2 ) keV component to have an emission 
measure which is < 30% of that of the 1.7 keV component. Employing eq.(jSJ), the volume 
filling factor of this component is then estimated as < [1 + (1.7/3.8) 2 /0.3] _1 = 0.6 (or 
equivalently rj c > 0.4). This is consistent with Figure El 

Figure [TO] also reconfirms that any cool emission below 0.7 keV is insignificant. The 
upper limit s are by more than a n order of magnitude below the prediction of the isobaric 
CF model (j Johnstone et al.lll992l ). if the model normalization is adjusted to the data at 1.7 
keV. Therefore, we conclude that the ICM around NGC 4696 is not cooling to temperatures 
much below ~ 0.7 keV. Analyzin g the RGS data from the present and an additional 110 ks 
observation, ISanders et al.l (120081 ) also confirmed that the emission measure from the central 
region decreases more steeply toward lower temperatures than in a simple CF picture. 



DISCUSSION 



5.1. Summary of the Obtained Results 



We analyzed the XMM-Newton EPIC data of the central region (< 12' or < 150 kpc) of 
the Centaurus cluster, assuming a spherical symmetry. The ICM has been characterized by 
three representative temperatures; Th ~ 3.8 keV, T c = 1.7 — 2.0 found within ~ 6' (75 kpc), 
and ~ 0.7 keV localized within ~ 1' (1 2 kpc) of the cen ter. The values of T c and Th are 
consistent with the ASCA measurements fllkebe et al.lll999l ). and sa tisfy eq.flTl). Similarly, th e 
detection of the 0.7 keV component reconfirms the Chandra result (ISanders &: Fabianll2002l ). 
The RGS data from the core region reinforced the presence of these discrete temperatures, 
and gave tight upper limits to emission with a temperature of < 0.5 keV. 



5.1.1. 2T vs. IT preference 

As far as the deprojected thin-shell spectra are individually analyzed, the data had no 
preference between the IP and 2P approaches. However, the two sets of EPIC spectra from 
3-dimensional thicker shells have revealed that the hot component is present even in a very 
inner (R < 5') region, and the cool component protrudes out beyond R ~ 3' ( §4.211 . The RGS 



result s from the cD region (§ 14. 4p provides an important support to this view. iTakahashi 



(120041 ) found that XMM-Newton spectra of the central regions of Abell 1795 also prefer the 



2P formalism to the IP modeling. 



In principle, the ICM could be in multi-temperature conditions (IKaastra et al.ll2004l ). 
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Since a spectrum resulting from such a su perposition of a range of temperatures is rather 
insensitive to the way of superposing them fjCraig fe Brown! 1 19761 ). an apparently 2P plasma 
may well be in a multi-temperature condition, and vise versa, particularly when the energy 
resolution is not too high. Indeed, the measured value of T c appears to distribute from 1.7 
to 2.0 keV, depending on the radius. Nevertheless, the RGS data preferred the presence 
of discrete temperatures, to a continuous temperature distribution ( §4.41) . Therefore, a 
multi-phase modeling of the ICM is considered inappropriate, at least the cD region of the 
Centaurus cluster is concerned. 

Based on these results, we suggest that the plasmas in the central ~ 75 kpc of the 
Centaurus cluster are better envisaged by invoking a few (two to three) discrete temperature 
components, rather than assuming a continuous radial temperature gradient or a condition 
of differential emission measure. These results support and extend our view of cD clusters 
developed with ASCA (Paper I). Below, we adopt the 2P view as our working hypothesis. 



When integrated up to R 



the 0.7 keV one have bolometric luminosities of L^° l 



12' (150 kpc), the hot component, the cool component, and 

3 x 10 43 erg s~\ L*° l = 1.0 x 10 43 erg 



s 1 , and Lq° 7 ' = 6 x 10 41 erg s 1 , respectively. According to llkebe et al.l (119991 ). L h 
by a factor of 2 if integrated further out to R = 30' (375 kpc). 



bol 



increases 



5.1.2. ICM abundances 



We have determined the radial abundance profiles of Fe, Si and O in the ICM, assuming 
that the different plasma components share the s ame abundances at a given radius. Those 
of the other elements are much more uncertain ( Takahashil 2004 ) . As fi rst discovered by 
Fukazawa et al.l (119941 ). and already reported by iMatsushita et al.l (j2007al ) by analyzing the 
present XMM-Newton data, the Fe and Si abundances have been reconfirmed to increase 
significantly toward the center (Figure [5]), and stay rather high even at the very center. 
The Si/Fe abu ndance ratio at the center, about 1.3 (Figure [5]), agrees with the ASCA 
measurements (IFukazawa et al.ll2000l ). 



In contrast to the behavior of Fe and Si, the O abundance is radially more constant, im 



plying a decrease in the O/Fe ratio toward the cen ter, in agreement with 



Matsushita et al. 



1998 



(|2007al). A s argued by various authors (Paper I ; iMushotzkv et al. 



2000: Sato et al. 2007a; Kawaharada 2006; Kawaharada et al. 



1996 



Fukazawa et al. 



20091 ) . this can natu- 



rally be understood as a result of varying contributions from type la and type II suparnovae: 
the outer-region ICM is contributed significantly by oxygen-rich products from type II su- 
pernovae, while that in the inner region by iron-rich products from type la supernovae that 
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occurred in the cD galaxy. 



So far, a fair number of other clusters have been known to exhibit abundance pro- 
files similar to those of the Centaurus cluster, with a central increase of iron and sil- 
icon, and a flat distrib ution of oxygen (e.g. 



2001 




Tamura et al. 


2001b 



Kawaharada et al.l 2009). These inc lude in partic ular recent Suzak u observations of the For 



nax 



cluster (IMatsushita et al. I l2007bh . Abell 1060 JSato et alil20~0~7bL and AWM 7 JSato et all 



20081 ). Therefore, the relative oxygen deficit in the cluter core region may be a common phe- 



nomenon. 



5.1.3. Evidence against Cooling Flows 



The present results argue against the CF scenario from several independent aspects. 
These include; the lack of cooling component ( §4.41) ; the lack of excess absorption ( §4. 3ft ; 
and the radial changes in the abundance ratios ( §5.1.2p which rule out large-scale inflows of 
the ICM. Below, we discuss the first two issues in some more details. 

The CF scenario relates the bolo metric lumino sity, Lcf, of a "cooling" portion of ICM 
with the mass deposition rate, M, as (lFabianlll994l ) 



5MkT h 
2/wrip 



(6) 



where fi is the mean molecular weight and m p is the proton mass. By substituting the 



measured value of L 



bol 



rbol 
"•^0.7 



1.1 x 10 43 erg s" 1 (§5.1) into L C f, we obtain M = 11 vr 



This in itself is consist ent with the previous results from ROSAT ( lEdge. Stewart Sz Fabian 
19921 ; iPeres et al.lll998l ) and ASCA (Paper I). However, we have found no evidence of emission 
from any plasma component with a temperature < 0.5 keV, with upper limits much tighter 
than the CF-predicted emission measure (Figure [TOl) . Therefore, the CF interpretation fails 
to explai n the strong cool em i ssion of the Centau r us cluster, as well a s in other clusters 
(Paper I; iTamura et al.ll2001al ; iKaastra et al.ll200ll ; iPeterson et al.ll200ll ). The value of M 
should not b e taken as an actual mass-deposition rate. 

We have shown in §4.31 that the previously reported excess absorption, at the center 
of the Centaurus cluster, is an artifact caused by the relative oxygen deficit. In fact, when 
the non-solar abundance ratios are properly considered, the absorption became consistent 
with being Galactic only, with any additional column at most several times 10 20 cm -2 . 
As discussed in §5.1.2[ the relative oxygen depletion in the cluster center appears to be a 
com mon feature. Then , the excess absorption, claimed previously in many other clusters 
(e.g. IWhite et al.l Il99ll ). will be explained away as well. Indeed, using the RGS spectra, 
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Peterson et al.l (120031 ) also confirmed the lack of excess absorption from several other clusters. 



Given these results, there is practically no evidence for the X-ray absorbing materials at least 
in nearby well-studied clusters. 



5.2. The 0.7 keV component 



Let us briefly consider the nature of the 0.7 keV component, which is localized to the 
cD galaxy, NGC 4696. Similar cool emission comonents are seen around other cD galaxies, 



including M< 


B7 in 


Fabian et al. 


2001a 



Churazov et al.ll200ll ; Molendi 1 12002 1 



Judging from the angular extent (~ 20 kpc) and the bolometric luminosity (6 x 10 41 erg 
s" 1 ), this 0.7 keV component could be hot interstellar medium (ISM) associated with the cD 
galaxy. Its temperature is consistent with the stellar velocity dispersion of NGC 4696, 274 
km s _1 . The plasma emitting this component (including bo th the central concentrat ion and 
the plume) has an estimated mass of a few times 10 9 M Q fjSanders Sz Fabian! |2002| ). which 
is also consistent with this ISM interpretation. Althoug h the estimated cooling time of this 
plasma is as short as ~ 10 8 yr (ISanders fc Fabian! |2002| ). the radiative energy loss could be 



supported by such processes in NGC 4696 as the supernova heating (a few times 10 41 erg 
s _1 ), and the accumulation of stellar mass loss which is dynamically "hot" on the galaxy 
scale. 

In contrast to the above arg uments, the complex filamentary morphology of this com- 



ponent (ISanders fc Fabian! |2002| ). including the plume like elongated tail structure, cannot 
easily be reconciled with the hot ISM interpretation, even considering interactions of the 
cD galaxy with the ICM. Sim i lar pr oblems with the cool gas in M87 (the Virgo cD galaxy) 
were pointed out by iMolendi I (120021 ) . who showed its close association with the synchrotron 
radio arms. Furthermore, as exploited in Paper I, the hot ISM interpretation could be ap- 
plied better to the cool (2 keV) component. Considering these, the 0.7 keV plasma could 
alternatively be associated with some magnetic structures around the cD galaxy, which in 
turn could be a result of past nuclear activity in NGC 4696. Further discussion of this issue 
is beyond the scope of the present paper. 



5.3. A Magnetosphere Model 



The successful 2P picture ( §4.21 §5.1.ip and the absence of significant CFs ( §5.1.31) in the 
Centaurus cluster suggest that its core region (within ~ 75 kpc) hosts a cool plasma phase 
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which is stably intermixed with the hot phase. While this reconfirms a prediction made in 
Paper I, it poses a series of questions to be answered; (i) how such a 2P configuration is 
realized; (ii) how the cool phase is heated against the radiative cooling; (iii) how the heating 
and cooling be balanced in a thermally stable manner; (iv) what determines the value of T c ; 
and (v) what produces the scaling of eq. ([!]). Below, we examine them one by one. 



5.3.1. Needs for ordered magnetic fields 



An immediate issue associated with the 2P scenario is (i) how to realize such a config- 
uration. We may suppose that the two phases are intermixed on a typical length scale of 
^ 10 kpc, as indicated by various sharp structures s een in X-ray images of some clusters 



(Fabian et al. 



2001b 



Machacek etHkood : Ilizukall2006h , including the "cold front" structure 



( IVikhlinin et al.ll200ll ). Then, employing the classical Spitzer heat conductivity, 



5 x 10~ 7 T 5/2 ergs s" 



cm 



-i K - 



(7) 



where T is the temperature in Kelvin, the time scale of conductive heat exchange between 
the two phases would be very short, as 



Tcond ~ 6 x 10" 



10 3 cm 3 



4 keV 



-5/2 



10 kpc 



yr 



( lSarazinlll988l ). The two phases w ould quickly become isothermal, unless the heat conduc- 
tion between them is suppressed f ITakahara fc Takaharalll979l ; lAsai. Fukuda fc Matsumoto 
2006h . 



The most n atural way to su ppress the heat conduction is to invoke magnetic fields, as 
argued, e.g., by iMolendi I (120021 ) to explain the cool plasma in M87. In fact, the ICM is 
known to be generally magnetized to a few /xG, and the aforementione d sharp X-ray feature s 
themselves are indicative of ordered strong (~ 10 /xG) magnetic fields (IVikhlinin et al.ll200ll ). 
In the particular case of the Cent aurus cluster, a m agnetic field strength up to ~ 25 /xG is 
reported from radio observations (ITaylor et al.ll2007l ). Although these field strengths are not 
necessarily higher than the equipartition value (~ 30 /xG), they are more than sufficient to 
completely suppress the heat transport across field lines, because the gyration radius of 4 
keV thermal electrons in a 1 /xG field, ~ 10 8 cm, is more than 10 orders of magnitude smaller 
than their classical mean free path for Coulomb scattering which also determines eq.([7]). In 
contrast, the heat conductivity along field lines may well be approximated by th e Spitzer 
value, eq.((7D, unless the field lines are highly tangled jMalvshkin fc Kulsrud IboOlh . 
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Based on solar analogy, ma gnetic field lines a round a cD galaxy may be classified grossly 
into "open" and "closed" ones (jMakishimal 1 199 71 ) . corresponding to coronal holes and coro- 
nae, respectively. This view relies only on a general and simple classification of the magnetic 
field-line topology, without invoking any particular magnetic configuration. Then, open-field 
regions, connected to the outer cluster volume, must be filled with the hot-phase ICM, and 
kept nearly isothermal by the efficient field-aligned heat conduction of eq.fjTj). In contrast, 
closed-field domains, connected to the cD galaxy and thermally insulated from the surround- 
ing hot phase, can take independent (possibly lower) temperatures of its own. In addition, 
the loop interior plasma is expected to be metal enriched by type la supernovae in the cD 
galaxy. 

From these considerations, we speculate that the plasma, filling numerous magnetic 
loops anchored to the cD galaxy, is observed as the cool phase of the Centaurus cluster (and of 
other similar objects). In other words, the cool phase ma y be regarded as a "magnetosphere" 
associated with the cD galaxy (Paper I; lMakishimalll997l ) . The magnetic fields, while working 
as a thermal insulator between the two phases, need not to be as strong as the equipartition 
value, since the loop-interior plasma is confined primarily by external pressure from the hot 
phase. This "magnetosphere" picture was already invoked successfully in Paper I to explain 
the ASCA results. 

Indeed, Ha-emitting filamentary features are optically observed in central r egions of 
some clusters on a spatial scale of a few to few tens kpc (e.g., IConselice et al.ll2Q0lh . and the 



brigh test of them can even be identified in some Chandra X-ra y images (e.g. 



Fabian et al. 



20031 ). The best example is the recent Chandra result on M87 flForman et al.ll2007t ). which 
reveals rich filamentary structures with a typical length of 10 to 50 kpc. In the Centaurus 
cluster as well, regions of strong magnetic fields are reported to be associated wi th filamentary 
structures that emit Ha photons and enhanced soft X-rays (ITaylor et al.ll2007l ). While these 
filam ents could be limb-b rightened edges of bubbles blown by the active galactic nucleus of 
M87 flForman et al.ll2007l ). they can alternatively trace magnetic structures such as assumed 
there. 

In a ddition to these X-ray and r a dio re sults, magnetohydrodynamic numerical simula- 
tions by lAsai. Fukuda fc Matsumotol (120071 ). including both the radiative cooling and the 
field-alighned heat conduction, show the emergence of low-temperature regions along the 
loop-shaped magnetic field lines. From these arguments, we believe we have found a viable 
answer to the issue (i) raised at the beginning of the present section. 
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5.3.2. The Rosner-Tucker-Viana mechanism 

We may assume that the cluster core region is fed with a constant heating lumiosity 
H, via some mechanism to be specified later in §5.41 In such a case, the volume heating 
rate will be generally proportional to the local ICM density, while the volume emissivity is 
obviously proportional to the density squared. Therefore, denser regions would preferentially 
cool and become even denser, leading to a thermal instability: this is the issue (iii). If, 
however, adopting the magnetosphere picture, the cool phase can be thermally stabilized 
by a built-in feedback mechanism , hereafter called RTV mechanism, originally developed by 



Rosner. Tucker fc Vaianal (119781 ) to explain quiet solar coronae. 



The RTV mechanism holds for a thin magnetic flux tube hydrostatically immersed in 
an external pressure, with the two loop footpoints anchored to a cool matter reservoir. The 
loop-interior plasma is assumed to be heated with a constant luminosity H, either uniformly 
along the loop or at the loop summit. The deposited heat H is assumed to flow toward 
the footpoints due to loop-aligned conduction of eq.([7j), and ultimately radiated away from 
various heights of the loop. Suppose that H was initially balanced by the radiative cooling, 
but some perturbation caused H to decrease. Then, the loop loses its internal pressure, and 
becomes thinner under the external compression, thus reducing the conductive heat flow 
along it. At the same time, some portion of the plasma flows into the matter reservoir, and 
the consequent decrease in the emission integral reduces the radiative energy loss. These 
concordant effects bring the loop into a new steady state corresponding to a lower value of 
H . Thus, the RTV mechanism thermally stabilizes the cool phase, and gives an answer to 
the issue (iii). 



Following the original work by lRosner. Tucker fc Vaianal (119781 ) , iKano fc Tsunetal (11995 



19961 ) used the Yohkoh spacecraft to examin the RTV mechanism. They confirmed that it 
is likely to be actually working in quiescent solar coronae, which keep emitting X-rays in 
an apparently steady manner for a time much longer than the nominal cooling time (~ 30 
minutes) . 



5.3.3. Cool-phase temperature 

The RTV mechanism further allows us to answer (iv), namely, to explain the absolute 
values of T c . Let us assume that the matter reservoir at the footpoints is much cooler than 
the loop interior, so that heat inflow to the reservoir can be neglected. Under such boundary 
conditions, the temperature of the loop-confined plasma is shown to take the maximum value 
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of 



T 

J- IT 



I.I - 1.4) x lO 3 ^/) 1 / 3 (K) 



(9) 



at the loop top (IRosner. Tucker fc Vaianalll978l ; iKano fc Tsunetalll995l ). Here, I (cm) is the 
loop half length, and p Gxt is the loop-confining external pressure in units of dyn cm" 2 . This 
i s so-called RTV tempe r ature scaling, which is actually confirmed in quiescent solar coronae 
f lKano fc Tsunetalll995l . Il996l ). The range of the coefficient, 1.1 to 1.4, reflects differences in 
the assumed heat input distribution along t he loop. The index of s caling in eq.Q depends 
weakly on the form of the cooling function (IKano fc Tsunetalll996l ): the particular value of 
1/3 is valid over a temperature range of 0.5-2 keV and near-solar abundances. 

Most importantly, T max in eq.Q is determined solely by the product p ex J, without 
depending on H. The loop-ineterior density n e becomes also independent of H, because 
p cxt oc n e T max is assumed to be constant. Changes in H affect only the loop cross section 
S, and hence the cool-phase X-ray luminosity which is by definition equal to H. More 
specifically, eq.flTJ) and eq.Q, together with H oc KT max S/l, can be combined into another 
scaling relation as 

SocHpZ'U-V* . (10) 

Thus, the loop becomes thicker under higher values of H. The loop also gets thicker if we 
decrease p cxt while keeping H constant, but in this case, T max automatically decreases via 
eq.®, so as to reduce the conductive heat flux (per unit cross section) oc KT max /l and hence 
keep H unchanged. 

In the present case of Centaurus, we may identify p ext with the hot-phase pr essure in the 
core region, namely p^ = In^kT^ ~ 1 x 10~ 10 dyn cm -2 with ~ 1 x 10~ 2 cm -3 (jlkebe et al. 
19991 ). and rewrite eq.Q as 



(2.0-2.5) 



Ph 



1 x 10~ 10 J V30kpc 



I 



1/3 



keV 



Here, we assumed that loop semi-lenght to be comparable to the narrower core radius, 
~ 30 kpc, of the d ouble-/3 modeling for the gravitational potential in the Centaurus cluster 
( llkebe et al.lll999l ). This is because each magnetic loop, with its interior having a higher 
density than the surrounding hot phase, may not become taller than the flat core radius of 
the central gravitational potential. The ass umed length is comparable to those of the soft 
X-ray filaments in M87 (IForman et al.l 120071 ). 



When such a loop is observed as an integrated entity, we expect to measure an X-ray 
temperature averaged over the loop length. Assuming the footpoint temperature to be 0.7 
keV ( §5.2|) . a simple numerical solution to the RTV equation predicts this average to be 
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~ 0.7 T max , or 1.4 — 1.8 keV from eq. ffTT]) . This is in an excellent agreement with T c actually 
measured in the Centaurus cluster, giving an answer to the issue (iv). 

In reality, the loop length may well scatter, with longer loops tending to be hotter 
according to eq. ffTT]) . Then, we expect the spherically averaged loop temperature to slightly 
increase outward, because each loop is hottest at its top region, and a larger radii would 
sample longer loops. This agrees with what we actually observe in Figure [5] 



5.3.4- The T c vs. T c scaling 

Taking it for granted that T c can be identified with ~ 0.7 T max , let us consider the issue 
(iv), namely, to explain the scaling of eq.fjlj). The luminosity of the hot phase in a cluster is 



empi rically known to scale with its temperature as oc ° (e.g.. IWhite. Jones fe Forman 
19971 ). On the other hand, we generally have oc n\ ^T h 1/2 , where R h is a typical radius of 



the hot phase. Elimination of from these two equations gives oc T^ 4 _R h 3 ^ 2 , and hence 
Ph oc n^Th oc T^ 4 R h 3 ^ 2 . Substituting this into eq. ffTT]) . we obtain 

T c oc T h 3/ X V \ /3 • (12) 

]/2 1/3 

Assuming that R h l h is relatively constant, T c indeed becomes roughly proportional to 
T h , providing an explanation to eq.fjTJ). 

As shown so far, the RTV mechanism can account for the stable coexistence of the cool 
and hot phases in the cluster core region. The scenario is independent of the amount of 
available heat input, because a lower value of H simply reduces the cool-phase luminosity. 
Furthermore, the RTV scaling can quantitatively explain the absolute values of T c , as well as 
its (near) proportionality to T h expressed by eq.fJT]). Thus, three of the four issues raised at 
the beginning of §5.31 have been answered. We therefore propose the magnetosphere model, 
incorporating the RTV mechanism, as a promising hypothesis to explain the plasma physics 
in the central regions of cD clusters in general. 



5.3.5. Galaxy vs. ICM interaction 



To complete our new viewpoint, we may ask ours elves why the bright c ool emission 
is observed predominantly around cD (or XD) galaxies (1 Jones fc Forman! 1 19841 . Figure 1 of 
Paper I), including NGC 4649. Conversly, non-cD galaxies, like NGC 4472 in the Virgo 
cluster, generally lack such bright cool X-ray emission, even if they are optically as luminous 
as their cD counterparts. 
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An outstanding characteristic of a cD galaxy is that it is nearly in a standstill at the 
bottom of the gravitational potential. Then, we may speculate that a substantial magne- 
tosphere can develop around a galaxy in a cluster environment only when it is at rest with 
respect to the ICM. When a galaxy is moving through the ICM, its magnetosphere will be 
disrupted by the ICM ram pressure, which can be one to two orders of magnitude higher 
than the magnetic pressure of the assumed magnetosphere (Paper I). Suc h strong interac 



tions are actually revealed by Chandra observations in nearby clusters (e.g.. iMachacek et al. 



20051 ; llizukal 120061 ). Furthermore, optical observati ons of distant clusters revea l numerous 



blue spiral galaxies with distorted morphology (e.g. 
suggest to be interacting with the ICM. 



van Dokkum fc Franxll200ll ). which we 



5.4. Possible Heating Mechanisms 



Having answered four out of the 5 questions raised at the beginning of §5.31 the only 
remaining issue is (ii), i.e., to seek for possible heating mechanisms to sustain the cool-phase 
luminosity. In the case of the Centaurus cluster, the target value is H = L^° l = 1.0 x 10 43 
erg s _1 . After Makishima (1997), Paper I, and the discussion conducted so far, we retain 
our magnetosphere viewpoint. 



The argument in §5.3.51 indicates that moving galaxies in a cluster interact with the 
ICM, and transfer their kinetic energies to the ICM. This process is estimated to proceed at 
a rate of 

L int ~ Nnm p v 3 7iR* nt (13) 



(ISarazinlll988l ). where N is the number of moving galaxies, n is the average ICM number 
density, v is the average galaxy velocity, and -R; nt is the interaction radius for each galaxy. 
We may use v = 1,010 k m s -1 which is a/3 times the vel ocity dispersion of the Centaurus 



cluster (586 km s 1 ) after iLucev. Currie fc Dickens I (119861 ). n ~ 0.8 x 10 3 cm 3 which is an 



average over the central 300 kpc (jlkebe et al.lll999l ). and N ~ 50. Furthermore, a relatively 
large value of R\ n t, say, ~ 5 kpc, is suggested by the argument in §5.3.51 an d by some X-ray 
observations (llizukal 120061 ). Through eq.t jTBl . these numbers yield L int ~ 4.5 x 10 43 erg s _1 , 
which well exceeds L^° l ~ lx 10 43 erg s _1 . Therefore, the galaxy-ICM interaction becomes a 
promising candidate for the ICM heating, and for the su ppression of CFs. These effects hay e 
been reproduced successfully by numerical simulations (lAsai. Fukuda fc Matsumotol 120071 ). 



How large is the kinetic energy stored in the moving galaxies? The stellar mass of 
the Centaurus clus ter within the cen tral 300 kpc, but excluding the cD galaxy, is estimated 
as 1.5 x 10 12 Mq (jlkebe et al.lll999l ) after correcting for the different values of the Hubble 
constant. We have then E^ n ~ 1.5 x 10 61 ergs, which is likely to be a lower limit because 
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the employed galaxy mass assumes a conservative mass-to-light ratio of 8. As a result, 
the time scale of the energy tansfer from the moving galaxies to the ICM is estimated as 
^kin/Lint ~ 1.1 x 10 10 yr. Therefore, the kinetic energy in the moving galaxies is large enough 
to sustain L int over the Hubble time. 

The galaxy-to-ICM energy transfer is expected to cause galaxies to gradually fall to the 
potential center. There are several pieces of observational evidence suggesting that this effect 
is actually taking place. In fact, the stellar-mass distributions in many clusters are known to 
be much more centrally peaked than that of the ICM. Furthermore, clear central decreases 
in "iron-mass to light ratio", observed from a fair number of clusters and groups (Paper I; 
Kawaharadal 120061 ; iKawaharada et al.ll2009l ). suggest that galaxies used to be distributed 
more widely than today, and polluted the ICM with metals out to the cluster edges. 

The energy released by the moving galaxies may be stored once in the I CM in the form o f 
turbulence and bulk flows. However, according to the Suzaku observation (jOta et al.l 120071 ). 
these kinetic energies in the ICM of the Centaurus cluster, as estimated from the width 
and position-dependent centroid shifts of the Fe-K line, cannot largely exceed its thermal 
energy. Therefore, the turbulent and bulk-flow energies must be dissipated and thermalized 
efficiently in the ICM, with a rate comparable to Lj nt . If a fair fraction of the dissipated 
energy is deposited onto the magnetosphere, the available heating luminosity of the cool 
phase becomes of the order of H ~ 1 x 10 43 erg s _1 which can just sustain the observed 
L^° l . We would rather say that L^° l is self-adjusted, via the RTV mechanism, to match 
the predetermined energy dissipation rate. The remaining portion of L int will heat the hot 
phase, to make it spatially more extended than the total gravitating mass. 

The actual dissipation of turbulence and bulk-flow energies may take place through 
magnetohydrodynamic processes. For example, moving galaxies will pick up magnetic field 
lines in the ICM and stretch them, exciting Alfven waves and causing field lines to reconnect. 
When the magnetic reconnection occurs between two closed magnetic loops, the released 
energy will be spent in heating the cool phase. That between closed and open field lines 
will expel a small potion of the metal-rich cool-phase plasma into the hot phase. The 
ejected plasmoid will become quickly isothermal with the hot phase, in terms of the electron 
temper ature via eq. flTD, as well as in terms of the ion ionization temperature (in a few times 
10 7 



yr; 



Masailll984T). The resultin g highly ionized heavy ions will remain there due to their 



slow diffusion (jEzawa et al.l 119971 ). and increase the metallicity of the hot phase. This can 
explain why the central metallicity increase is observed in both phases. 

In addition to the above mechanism, there can be another ICM heating mechanism 
localized around each cD galaxy; namely, a shrink of its self-gravitating core. Since this has 
already been discussed in Paper I, we simply mention that this mechanism can also account 
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for a heating luminosity of ~ 1 x 10 43 erg s 



6. CONCLUSION 

Analyzing the XMM-Newton EPIC and RGS data of the Centaurus cluster, we have 
obtained the following results. 

1. The ICM in the central ~ 75 kpc can be described better by the two-phase view (plus 
the 0.7 keV component at < 12 kpc) than by the single-phase picture. 

2. There is no evidence of cooling flows. 

3. The iron and silicon abundances increases significantly in the central region, while that 
of oxygen is radially constant. 

4. The previously reported excess X-ray absorption disappears when the central oxygen 
deficit is considered. 



5. The overall results agree with those by llkebe et al.l fll999l ) who used ASCA and ROSAT, 



To explain these observations, we have developed a working hypothesis (after Paper I), 
in which the cool phase is thought to constitute a magnetosphere associated with the cD 
galaxy. Incorporating the Rosner- Tucker- Vaiana mechanism, it can account for the following 
essential features of the Centaurus cluster, as well as of similar cD clusters. 

1. The hot and cool ICM phases co-exsit in the cluster core region. 

2. The ICM is provided with a heating luminosity, which is high enough to sustain the 
X-ray emission against the radiative cooling. 

3. The cool phase is kept thermally stable. 

4. The cool phase has a typical temperature of ~ 2 keV, with a mild outward increase. 

5. A good proportionality holds between the cool and hot temperatures. 

6. The ICM metallicity becomes enhanced in both phases toward the center. 

The present results thus provide a new insight into the physics of clusters of galaxies. 
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Apendix A: The Deprojection Matirx 

This Appendix is meant to provide numerical details of the deprojection process de- 
scribed in §3.2.31 Let us consider 11 annular regions, of which the outer radii are at 0'.5, 
l'.O, 1'.5, 2', 3', 4', 5', 6', 8', 10', and 12', and the corresponding 11 shell regions. These 
divisions are the same as used in Figures [2J HJ and except that the innermost 4 regions 
therein are here re-arranged into two regions just for clearer presentation. The reconstructed 
"Shell C" spectrum, namely C = £j =3 Sj, is not affected by this simplified treatment of inner 
annuli/shells. Let (A±, A%, An) be the projected spectra from the 11 annular regions at a 
given energy, while (Si, S 2 , Sn) those from the 11 shells. Note An has the outer boundary 
at 12', while Sn includes emission outside that. Then, let Aiij denote the projection matrix 
which converts {Si} into {Aj}. In the present geometrical setting, it is numerically given as 



/ A, \ 

A 2 
A, 
Ai 
A 5 
A e 
A 7 
A s 
A, 
A 10 
V A u 



J 



( 1.000 










V o 



0.258 
0.742 





















0.083 
0.329 
0.588 



















0.041 
0.137 
0.321 
0.501 

















0.020 
0.063 
0.116 
0.213 
0.588 















0.010 
0.031 
0.055 
0.083 
0.321 
0.501 













0.006 
0.019 
0.032 
0.047 
0.151 
0.303 
0.443 











0.004 
0.013 
0.021 
0.031 
0.093 
0.152 
0.285 
0.401 









0.003 
0.008 
0.013 
0.018 
0.055 
0.083 
0.122 
0.199 
0.501 







0.002 
0.005 
0.008 
0.011 
0.032 
0.047 
0.064 
0.086 
0.303 
0.443 




0.004 \ 




( St 


0.011 




s 2 


0.018 




s 3 


0.026 






0.074 




s 5 


0.105 






0.137 




S7 


0.174 




Ss 


0.477 




Sq 


0.724 






1.113 ) 







On the assumption that the emission is uniform within each individual shell, these matrix 
elements are uniquely determined by the geometry alone, without any dependence on the 
emission model or instrumental effects. The only exception is the last column, namely Aii t u 
{i = 1,2, ..,11), which assumes the emission outside 12' as explained in §3.2.31 Although 
the matrix elements are shown with only 3 digits below decimal points, this is again for 
simplicity; the calculation is preformed with double precision. 

By inverting this matrix, we can calculate the deprojection matrix {T>ij} defined in 
eq.©. Numerically, it is calculated as 
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From eq. (T151) , the Shell C spectrum C at each energy can be expressed as 

C = (1.700A 3 +0.909A 4 +1.036A5+0.999A 6 +1.003A 7 )-(1.493A 8 +0.005A 9 +0.085Aio+0.045A 11 ). 

(16) 

Thus, the outermost annulus An, which is most subject to the background uncertainty as 
well as to that of the emission outside R = 12', contributes geometrically only less than 5% to 
the Shell C spectrum, compared to the individual shells A 3 through A-j. In addition, the flux 
in An is ~ 3 to ~ 30 times lower than those in the 3rd to 7th annuli (Figure [2]). Considering 
these, we estimate that the emission from the R > 12' region contributes physically by no 
more than ~ 1% to the Shell C spectrum. Therefore, any systematic error associated with 
An is concluded to be negligible ( §4.2.61) . 

Another exercise to be conducted using eq. f|T5|) is the issue raised in § 14.2.11 i.e., 
how the assumption of emission constancy within each thin shell affects our Shell C spec- 
trum. As mentionned there, let us examine what happens if we purposely neglected spec- 
tral differences between A 5 and A 6 (the 3rd and 4th thin shell constructing Shell C, re- 
spectively), where the IP properties change rather steeply with R. Taking the Fe-K line 
equivalent width as an example, we may write the fluxes in the Fe-K line energy region as 
A 5 = B 5 (l + W 5 ) and A e = B e (l + W e ), where B represents continuum and W the line equiv- 
alent width. Using the difference 5W = W& — Wq and the emission-measure-weighted mean 
W = {B 5 W 5 + B 6 W 6 )/(B 5 + W 6 ), we may rewrite as A 5 = B 5 [l + W+ B (i 5W/(B 5 + B 6 )] 
and A Q = B Q [l + W - B 5 5W/(B 5 + Bq)] . Substitution of this into eq.f^B) readily yields 

C = C + (1.0365s + 0.999fi 6 )(l + W) + 0.037a (B 5 + B 6 )5W (17) 

where C represents terms unrelated to A 5 or Aq, while a = B 5 Bq/(B$ + Bq) 2 < 1/4 is a 
numerical factor. Thus, the negligence of SW between the two annular spectra has only an 
effect of at most ~ 0.0151U in the Shell C spectrum. The essence is that the deprojection 
procedure has mathematically "differential" nature exaggerating small-scale features, while 
this problem is largely removed by constructing a thick shell wherein a process of integration 
is involved. This is evidenced by eq. ffTBT) . where C is close to a simple summation of ^3 
through A 7 , except that A 3 must be given a 1.7 times larger weight and A s (the shell just 
outside Shell C) must be subtracted with a relatively large weight of 1.5. 

Apendix B: A Numerical Simulation of an IP Cluster 

In order to assess the validity of our "thick shell" method, the following numerical 
simulation was performed. We started from the best-fit IP numerical models (Figure H]) 
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obtained for the 13 deprojected thin shell regions which are defined in Figure [2j Then, 
to represent the radial spectral changes more smoothly, each of the 5 shells constituting 
Shell C (covering r = V — 5') was subdivided into inner and outer halves, and appropriate 
temperature and emission measure were assigned to the inner/outer pair. These model 
spectra, now comprising 18 shells, were projected numerically onto the corresponding 18 
annuli, to form numerical spectral models for the 18 annular regions. At this stage, the 
inner /outer pair were numerically recombined together. The derived 13 model spectra were 
convolved with the EPIC energy responses, and were given Poisson errors simulating the 
actual observation. Then, in the same manner as the actual data, we analyzed these 13 
annular spectra via deprojection. 

Figure [TT] shows the simulated Shell C spectra, compared with the "sysnthetic IP spec- 
tra" introduced in §4.2.21 The comparison is fully satisfactory with x 2 jv = 324.96/322 = 
1.009. For reference, a 2P fit with two APEC components, in the same manner as explained 
in §4.2.31 gave x 2 / u = 312.9/309 = 1.013, which is similarly acceptable. In other words, 
the IP and 2P interpretations degenerate unlike the case of the actual data. Furthermore, 
unlike the 2P fit to the actual Shell C spectra which revealed the Th ~ 4 keV component, 
the 2P fit to the simulated data yielded T c = 2.00+°;^ keV and T h = 3.11+°;^ keV, which 
are both consistent with the assumed IP temperature range (2.03 to 3.30 keV). 
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Fig. 1. — EPIC PN (panel a) and M0S1 (panel b) light curves in 2-7 keV from the present 
observation, obtained by excluding the central (r < 8') region. The bin width is 100 sec. 
The dashed lines indicate ±2<r levels from the mean values in the quiescent periods. 
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Fig. 2. — Background-subtracted EPIC-PN spectra of the Centaurus cluster extracted from 
annular regions (panel a), and their ratios to that averaged over the region of 5' < r < 12' 
(panel b). Outer radii used to extract the spectra are [a] 10", [b]25", [c]40", [d]l', [e]l'.5, [f]2', 
[g]3', [h]4', [i]5', [j]6', [k]8', [1] 10', and [m]12'. The spectra are not corrected for the vignetting. 
In panel (a), only 5 out of the 13 spectra are shown for clarity. The ratios in panel (b) are 
vertically offset to avoid overlap. 
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Fig. 3. — Deprojected PN spectra (colors) from concentric thin shells, compared with their 
non-deprojected counterparts (black). The regions used are; [cd]0'.5 — 1', [f]l'.5 — 2', [i]4' — 5', 
and [k]6' — 8'. To avoid heavy overlaps, a series of spectral pairs are shown on every other 
panel. 
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Fig. 4. — Radial profiles of (a) the temperature, (b) reduced chi-squared, and the abundances 
of (c) Fe, (d) Si, and (e) of the Centaurus cluster, obtained with the IP modeling. A 
~ 0.7 keV component is added within r < 1'. Filled circles and open diamonds represent 
results from the non-deprojected and deprojected spectra, respectively. 
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Fig. 5. — The same as Figure HI but using the 2P model. The results outside 6', where 
the cool component is not required, are th e same as those o f Figure HI In panel (c), the 
Fe abundance profile obtained with ASCA (jlkebe et al.lll999l ) is shown in the dashed line, 
together with its typical error size (30 %). 
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Fig. 6. — (a) Radial profiles of the emission measure per unit volume of the hot (red dia- 
monds), cool (blue crosses), and the 0.7 keV (green squares) components, derived through 
the 2T fits to the deprojected EPIC spectra. Dashed and dotted horizontal lines indi- 
cate the radial ranges covered by Shell C and Shell P, respectively, (b) Radial profile of 
the vol ume filling factor of the cool component, calculated via eq.flHj). That obtained with 
ASCA fllkebe et al.lll999l ) is also shown with a dashed line, together with its typical error 
(25%). 
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Fig. 7. — (a) Actual spectra from the thick Shell C, compared with the synthetic IP model 
spectra (solid histograms) constructed from the best-fit IP models to the constituent 5 thin- 
shell spectra (see text §4.2.21) . The PN, MOS1, and MOS2 spectra are presented in black, 
red, and green, respectively, (b) The same Shell C spectra, fitted jointly with a simple 2P 
model of free temperatures and free abundances ( §4.2.31) . The fit goodness is summarized in 
Table 1. 
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Fig. 8. — A temperature map of the Centaurus cluster, derived by analyzing the projected 
EPIC spectra in four azimuthal sectors and the 12 radial intervales. The MOS and PN 
spectra are fitted simultaneously by a IP model (§ 14.1.21) . The map scale is in units of 
arcsec, with north up and east to the left. The temperature color scale is in the unit of keV. 
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(a) RGS1+RGS2, 1st + 2nd orders 
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Fig. 9. — (a) Background-subtracted RGS flux spectrum made by combining four spectra, 
namely, the 1st and 2nd order spectra from the two RGS units. The positions of the H-like 
and He-like K-lines of several major elements are indicated with dashed and dotted lines, 
respectively. Color tick marks indicate Fe-L lines of various ionization states, (b) RGS1 
spectra, fitted jointly with a 2P model. The first order spectrum and its model prediction 
are shown in black and red respectively, whereas the 2nd order data and its model predictions 
are in green and orange, respectively, (c) The same as panel (b), but for RGS2. 



-50- 




0.3 0.5 0.7 1.0 2.0 3.0 4.0 

Plasma temperature (keV) 



Fig. 10. — Distribution of the normalizations of individual temperature components, ob- 
tained with the multi temperature fit to the RGS data. The curve in dicates the prediction b y 
the isobaric cooling flow model, for plasmas of 1 solar abundances (I Johnstone et al.lll992l ). 
The vertical scale is arbitrary. 
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Fig. 11. — Spectra from Shell C of a simulated IP cluster, compared with the convolved 
numerical models that are used to define the 5 constituent annular spectra (indicated only 
for MOS1 with different colors). The PN, MOS1, and MOS2 spectra are presented in black, 
red, and green, respectively. Their appearance differs from that of Fig. simply because 
the energy- dependent effective areas, which are once removed in deprojection, are restored 
in Fig. [7] while not here for simplicity. 
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Table 1: Systematic uncertainties (1 a) assigned to the EPIC background data. 





band (keV) 


errors 


note 


PN 


< 1.35 


8% 


diffuse X-rays 




1.35-1.6 


8% 


Al line 




1.6-7.3 


3% 






7.3-9.2 


8% 


Ni, Cu, and Zn lines 




> 9.2 


3% 




MOS 


< 1.35 


8% 


diffuse X-rays 




1.38-1.85 


8% 


Al and Si lines 




> 1.85 


3% 





Table 2: Summary of the fit goodness to the Shell C spectra. 



Model 


Condition 


2 

X 


V 


X 2 /" 




IP 


synthetic, no adjust. 


739 


378 


1.96 


4.2.2 




norm, adjustted 


724 


373 


1.94 


4.2.2 




revised sys. error 


418 


378 


1.11 


4.2.4 


2P 


simple 2P fit 


688 


363 


1.90 


4.2.3 




synthetic, no adjust. 


678 


378 


1.79 


4.2.3 




revised sys. error 


381 


378 


1.01 


4.2.4 
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Table 3: Results of the model fit to the RGS spectra. 



Parameter Comp. 1 Comp. 2 
Temperature 1.70+^ 0.77l^ 

Gaussian a Q5".8%j Ig^ggj 
Abundances (solar) 

O 0.39 ±0.04 

Ne 0.60 ±0.12 

Mg 0.75±g;il 

si i.09±g"il 

Fe 0.75j™ 

Ni LOl^j 

X 2 /d.o.f. 985/727 



